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Abstract 

The perimeter and area generating functions of exactly solvable poly- 
gon models satisfy (/-functional equations, where q is the area variable. 
The behaviour in the vicinity of the point where the perimeter generat- 
ing function diverges can often be described by a scaling function. We 
develop the method of q-linear approximants in order to extract the ap- 
proximate scaling behaviour of polygon models when an exact solution 
is not known. We test the validity of our method by approximating ex- 
actly solvable q-linear polygon models. This leads to scaling functions 
for a number of q-linear polygon models, notably generalized rectangles, 
Ferrers diagrams, and stacks. 

1 Introduction 

Models of polygons and related combinatorial objects have received considerable 
attention in recent years (for a recent monograph, see j8|). They are of interest 
in physics as models of vesicles or polymer molecules in solution. The interplay 
between bulk energy and surface energy in these models gives rise to a phase 
transition from an extended phase to a compact, ball-shaped phase |^|. There 
have been many studies of combinatorial aspects of these models, including 
a general method for deriving the perimeter and area generating function of 
column-convex models S. Less is known about analytic aspects of the solutions, 
which are needed to understand the phase transitions of these models. Scaling 
functions which describe the crossover behaviour at critical points have been 
computed for a number of polygon models, mostly by indirect methods such as 
from a semi-continuous version of the models |l4|, . The only direct derivation 
of scaling behaviour has been for staircase polygons |H| by methods of uniform 
asymptotic expansions. There is, however, no known general method to obtain 
scaling functions directly from functional equations. 

This paper presents such a method in the simplest case of a g-linear func- 
tional equation in the perimeter variable. This class of functional equations is 
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satisfied by rectangles, Ferrers diagrams, and stacks Q. As a step towards the 
analysis of more complicated classes, we introduce g-linear approximants of first 
order so as to analyze models which do not obey a first-order g-linear equation, 
but which can be well approximated by one. We will test our method by ap- 
proximating exactly solvable g-linear polygon models of generalized rectangles, 
Ferrers diagrams and stacks. In particular, we will analyze the model of Ferrers 
diagrams with a hole and obtain a differential equation for the scaling func- 
tion by analysis of the approximants. We discuss the connection between this 
new type of approximant and the method of partial differential approximants 
[Q [It], [l^, |l8). Finally we indicate how our methods can be extended to more 
general classes of polygon models. 

In a subsequent publication we will consider g-quadratic and other non-linear 
approximants, which we expect will give good approximations to the scaling 
function of as yet unsolved models, such as self-avoiding polygons. 

2 Phase diagrams and scaling functions 

Let us briefly review phase diagrams of g-linear polygon models in order to 
fix our notation. (We follow JlJ, |ll| [ll]]). The perimeter and area generating 
function of a polygon model is given by 

oc oo 

f{x,y,q)= ^2 fr,s,nx r y s q n = ^2fn(x,y)q n , (1) 

r.s,n— 1 n— 1 

where f r , s ,n denotes the number of configurations of area n, horizontal perime- 
ter r and vertical perimeter s. We introduce the area activity q, the horizontal 
perimeter activity x, and the vertical perimeter activity y. The perimeter gen- 
crating function of the polygon model is given by f(x, y, 1). A phase diagram 
is the graph of the radius of convergence of f(x,y,q) in the parameter space 
x, y, q. Let us consider the isotropic version f(t, q) :— fit, t, q) of the model, with 
t denoting the total perimeter, so that the phase diagram is two-dimensional. 
For a typical g-linear polygon model such as Ferrers diagrams or stacks, defined 
below, the phase diagram is as depicted in Figure [|. Let us interpret the 
phase diagram using the grand-canonical ensemble in which we count, for fixed 
area, all polygons by perimeter. The curve q c (t) where the polygon generating 
function diverges is related to the free energy per unit area of the ensemble, 

- log q c {t) = lim -log/ n (t). (2) 

n — >oo ti 

The phase q c (t) < 1 consists of inflated polygons, whose perimeter grows like 
their area. The phase q c {t) = 1 consists of ball-shaped polygons, their perimeter 
growing as the square-root of their area. This results in a vanishing free energy 
for the ensemble. This behaviour is characteristic of a first order phase transition 
at the point tf where both phases meet. In the ball-shaped phase q c {t) = 1, there 
are contributions to the boundary free energy, howeverQ Let us denote by t c the 

1 The boundary free energy is defined as the limit fb(t) = lim n _,oo -h= log f n (t). 
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Figure 1 : Phase diagram of a typical q- linear polygon model 



point where the perimeter generating function diverges. At this point a phase 
transition in the boundary free energy occurs: For t < t c , the contributions to 
the boundary free energy are given by polygons of finite size, whereas for t > t c 
the contributions to the boundary free energy derive from polygons of infinite 
size. 

In the remainder of this paper we will concentrate on the critical behaviour 
about the point t c where the perimeter generating function diverges. This point 
is the natural one to look at from the perspective of power series approximations, 
which we employ. Moreover, for rectangles and more complicated models such 
as self-avoiding polygons, the distinction between the two phase transitions is 
irrelevant since t c and tf coincide for these models. 

To describe the singular behaviour about t c in more detail, consider f(t, q) 
for t fixed, as q approaches unity. For (/-linear polygon models, q — 1 is a point 
of an essential singularity in the generating function: For t < t c , f converges to 
a finite limit. If t = t c , f has a power-law divergence with an exponent generally 
different from that of the perimeter generating function. If t > t c , f diverges 
with an essential singularity. In many cases, the crossover between these types 
of critical behaviour can be described by a scaling function P(s) of combined 
argument s — (t c — t)(l — g) _< ^, 

The asymptotic behaviour of the scaling function at infinity is related to the 
behaviour of f(t,q) for t < t c . To see this, assume that f(t,q) admits an 
asymptotic expansion of the form 

oo 

f{t,q) = Y J fn{t)(l-q) n (t<t c ) (4) 

about q = 1, where the leading contributions of the coefficients f n {t) are given 
by 

fn(t) = u \, + o ^- t y n+1 )^ ( 5 ) 



3 



as t approaches t c ^ For g-linear polygon models, the asymptotic expansion 
can be computed recursively from the defining functional equation, ft can be 
inferred from (0) that the existence of a scaling function implies the restriction 

On 

on the exponents j n . Moreover, it can be seen that the numbers p n are the 
coefficients in the asymptotic expansion of the scaling function 

oo 

P(s) = X>n*r 7 ". (7) 

n=0 

We assume that the scaling function P(s) is regular at the origin. (This as- 
sumption is not always fulfilled. The simplest counterexample is the model of 
rectangles in its isotropic version.) In this situation, the behaviour of / at t = t c 
is given by 

(rr|« fa- 1- ). (8) 

The exponents 9 and 70 are called critical exponents of the model. 6 describes the 
behaviour of f(t c , q) about q = l, whereas 70 describes the power- law behaviour 
of the perimeter generating function about t c . The exponent <j) is called the 
crossover exponent and relates the two critical exponents, see (|^). 



3 (/-linear polygon models 

We call a polygon model g-linear of iVth order if its generating function satisfies 
a q-linear functional equation^] 

N 

G{t, q)=J2 a ^> l) G (l kt > l) + 6 (*> ?). ( 9 ) 
fe=l 

where ak(t, q) and b(t, q) are rational functions in t and q. Here, t may denote the 
total or horizontal or vertical perimeter. Explicit realizations include rectangles, 
Ferrers diagrams, and stacks flli|j . They all satisfy a g-linear functional equation 
of first order in the horizontal perimeter activity x. 

We may construct new polygon models from given ones by allowing for 
decorations of the polygons. In this way we may obtain models of polygons 
with holes, of coloured polygons, or the like. Decorations may be interpreted 
in physical terms as allowing for a refined structure of polygons, which may be 
a better approximation to vesicles than undecorated polygons. These models 
have a natural interpretation in terms of random tilings of polygonal shape ||l6f . 

2 The coefficients may have a different asymptotic form, see rectangles at y = 1 below. 
3 Properties of solutions of g-linear functional equations if q < 1 have been studied in [El, 
for recent results see lp|. 
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The question arises as to which decorations lead to models which continue 
to satisfy q-linear functional equations. This is the case for models whose gen- 
erating function can be obtained by application of a linear differential operator 
(w.r.t. x, y, and q) to the generating function of the undecorated model. For 
example, the model of Ferrers diagrams with a Thole, defined in the appendix, 
satisfies a q- linear functional equation of order 3. Let us now focus on the 
simplest class of decorated polygon models which satisfy a g-linear functional 
equation of first order. To this end, consider first decorated rectangles of unit 
height. We denote the generating function of this model by b\(x, y, q). For ex- 
ample, the generating function of a rectangle with exactly k black unit squares 
is given by 

y(qx) k 



bi(x,y,q) 



(10) 



while the generating function of undecorated rectangles is given by bi(x, y, q) = 
yqx/(l — qx). Let us now consider the polygon model of rectangles with a 
decorated top layer. As indicated in Fig. 0, we can consider all rectangles of 
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Figure 2: Graphical representation of the functional equations for rectangles, Ferrers 
diagrams, and stacks with a decorated top layer. 



height to + 1 as being obtained from rectangles of height m by adding a row 
of squares below the bottom layer. This construction misses out all decorated 
rectangles of height 1. A similar construction can be applied to Ferrers diagrams 
and stacks with a decorated top layer. The figure indicates that the models 
defined above satisfy a g-linear functional equation of first order 

y 

G s {x,y,q) = — — — G s (qx,y,q) +bt(z,y,q), (11) 
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where s = 0, 1,2 denotes rectangles, Ferrers diagrams, and stacks respectively, 
and b%(x, y, q) is the generating function of decorated rectangles of unit height. 
Equation ( |TT| ) can be iterated to give a closed form for the area and perimeter 
generating function 

Gs(x ' y ' 5)= S ' (12) 

where (t;q) n — IIfc = o(l~ q k t) denotes the g-product. The perimeter generating 
function can be obtained from the functional equation by setting q = 1 and 
solving for G s (x,y, 1), 

G s (x, y, 1) = n (1 ~f S h(x, y, 1). (13) 
(1 - x) s - y 

The models of rectangles, Ferrers diagrams and stacks are recovered as special 
cases where the trivial decoration is chosen. 



4 Scaling functions via dominant balance 

We now discuss how to derive scaling exponents and scaling functions from q- 
linear functional equations about the critical point where the perimeter generat- 
ing function diverges. Our technique relies on the method of dominant balance 
[pi , which has previously been applied to the derivation of scaling equations for 
a number of semi-continuous models 13, 14|]. Our approach lies in deriving scal- 
ing functions for the discrete models by manipulating the g-functional equation 
directly. It consists of three steps: Firstly, the critical point is shifted to the 
origin by a change of variables. Then a scaling variable is introduced, and a 
consistent set of scaling exponents is sought. Finally, the resulting differential 
or difference equation is solved. 

Consider the q-linear functional equation of first order 

o-o{x, q)f{x, q) - ai(x, q)f(qx, q) - b(x, q) = 0. (14) 

Bearing later applications in mind, we restrict ao(x,q) 1 ai(x,q),b(x,q) to poly- 
nomials in x and q. For readability, we suppress all subsequent dependencies 
on q. As a first step, assume that the critical point is at x — x c and q = 1. 
We expand the functional equation about the critical point. To this end we 
introduce small variables e = 1 — q and s — x c — x, and define P(s) = f(x c — s). 
In these variables, the functional equation reads 

(oo - oi) (x c -s)P(s)- ai (x c -s) J ^r p ( s )-^~ s ) = °- ( 15 ) 

n—1 



Second, we introduce scaled quantities 



and write equation (|15|) in terms of s and P. These quantities are just the 
scaling variable of combined argument and the scaling function (|^) . 

Here, <fi is assumed to be positive. This leads to additional factors of e in 
each summand. The scaling equation results from taking only the terms with 
smallest exponents in e, for suitable choices of the exponents and 4>. 

Let us analyze the contribution from the n-th order in the expansion of f(qx). 
Scaling leads to exponents of the form n(l — <f>). This implies the constraint 
(f> < 1, since other values lead to arbitrarily small exponents with increasing 
n. If <f> < 1, only the first order contributes, leading to x c j^P(s). If <fi = 1, 
higher orders cannot be ignored, and the dominant part of the sum equals 
P(s + x c ) — P{s). This results in a differential or difference equation for the 
scaling function. We will present typical examples of both kinds below. 

We concentrate on critical points given by the smallest pole of the perimeter 
generating function, that is at the smallest positive x c satisfying 

ai(x c , 1) = a (x c , 1). (17) 

For polygon models, the coefficients of the perimeter generating function are all 
positive. This generating function is obtained by setting q = 1 in ([l5]), and hence 
implies that b(x) / (ao(x) — ai(x)) has nonnegative Taylor-coefficients. This im- 
plies in particular that b(x c , 1) is non-zero. Therefore, the leading contribution 
of b(x) in the scaling limit is of order e°. The contributions from ao(x) — ai(x) 
and from a\{x) have to be analyzed for each model separately. To obtain a 
nontrivial scaling equation, we demand that we get contributions from each of 
the three terms in (|l5|). This means that all three exponents in e have to be 
equal. We thus arrive at a set of equations determining 9 and </>, 

...-6 = ...-9+{l-<i)) = 0, (18) 

where . . . denotes contributions from ao(x) — ai(x) and ai(x), respectively. We 
discuss particular examples below. We finally mention that a scaling analysis 
of the g-linear equation of order N can be carried out by the same method. If 
4> < 1, this results in a differential equation, as found above for the case N = 1. 

4.1 Ferrers diagrams 

Ferrers diagrams satisfy a q-linear functional equation with 

ao(x,q) = 1-qx, (19) 
ai(x,q) = y, 

K x ii) = yi x - 

The perimeter generating function diverges at x c — 1 — y. ai(x,q) and b(x,q) 
are nonzero at the critical point, whereas ao(x, q) — cii(x, q) vanishes linearly in 
s — x c — x. The leading contributions in the three terms of equation (|15|) give 



■-0 = -0+(l-0)=O. (20) 
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This leads to exponents 

e = \, ^4 (2i) 

and to the differential equation 

1 -sP{s) = ^P(s) + 1. (22) 



y{\ - y) ds 

The solution, which is uniquely determined if we demand power-law behaviour 
as s — > oo, is given by a complementary error- function 

lira „ / s \ ( s 2 \ 



where a = y(l — y). The asymptotic behaviour is given by 



P(s) ~ - (s^+oo) (24) 
P(s) ~ y/2na exp ^— ^ (s — > — oo) 

P(0) 

This is in agreement with both the scaling function found for the semi-continuous 
model and with the asymptotic behaviour computed previously jL4] . 

4.2 Decorated rectangles 

We next consider the model of rectangles with a decorated top layer. The 
decoration consists of exactly k black squares. The generating function of the 
top layer is 

y(qx) k 

(1 - qxf 

As shown above, the model satisfies a (/-linear functional equation with 



h(x,y,q)= (25) 



a (x,q) = (l-qx) k+1 , (26) 
ai(x,q) = y{l-qx) k+1 , 
b(x,q) = y{qx) k . 

The perimeter generating function has a pole of order k + 1 at x c = 1. In 
contrast with the model of Ferrers diagrams, this point coincides with the point 
Xf where the first order phase transition in the free energy occurs. The case 
k = is closely related to rectangles^: Comparison of (j|^) with the functional 
equation for rectangles shows that both models obey the same equations about 
the critical point and hence have the same scaling functions. 

4 For rectangles, bi(x,y,q) = yqx/(l — qx). 
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The leading terms in the three summands of equation ( |15[ ) give 

(fc + l)0-0 = (fc + l)0-0+(l-0)=Q. (27) 
This leads to exponents 

9 = k+l, <j)=l, (28) 
and we get the difference equation 

(1 + s) k+1 (P k (s) - yP k (s + 1)) - y = 0. (29) 

In order to obtain the correct scaling function from this recursion, we have to 
force the asymptotic behaviour as s — > oo to be of power-law type. This fixes 
the constant term, and for y < 1 we arrive at the Lerch functions 

oo „ 

The asymptotic behaviour is given by 

Pk(s) ~ T^i^TT (s^+oo) (31) 



Pk(0) = £ 



y s n 

y r, 



n 

n=l 



fc+1 ' 



The case k = (rectangles) is in agreement with both the scaling function of 
the semi-continuous model and the asymptotics computed in Oj. 

The case y = 1 is different. The perimeter generating function diverges for 
all values of x, but it is possible to obtain scaling behaviour about x c = 1, as q 
approaches unity. For k > 0, the solution of the difference equation is given by 

(-l) k+1 

Pk ^ = kl + (32) 

where vE , fc(x) denotes the fc-th derivative of the ^-function ^(x) = d x logT(a;). 
The asymptotic behaviour is given by 

P k (s) ~ ~ (s^+oo) (33) 

P k (o) - C(* + i), 

where C(fc) = J2^Li n ~ k denotes the Riemann zeta function. 

If k = 0, the scaling function diverges logarithmically as s — > oo. In order 
to obtain a well-defined crossover behaviour, we compensate for this by the 
addition of a logarithmic term, 

P (s) = -*(l + s)-loge. (34) 
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The asymptotic behaviour then follows as 



Ms) 
Po(0) 



log(se) = 
loge + 7, 



logs 



(s — ► +00) 



(35) 



where 7 = — J °° e - ' lni alt « 0.57721 denotes Euler's constant. This gives the 
(non-uniform) asymptotic behaviour first computed in Jl4| . The scaling function 
describing this behaviour was not previously known, however. 

5 g-linear approximants 

We now consider the more general situation where we cannot obtain the gener- 
ating function, but only a finite number of terms thereof. 

5.1 The method 

The basic idea of q-linear approximants (of first order) is to fit a q-linear func- 
tional equation to a function /(t, q), 



such that (|3^) is exact up to a given degree in the variables. Here, we restrict 
ao{t,q), ai(t,q), and b(t,q) to be polynomials in t and q of degree n t and n q , 
say. The coefficients of the polynomials can be found by solving the system of 
linear equations deriving from the expansion of ( p6| ) in its two variables. This 
process is not unique, since there are many choices of sets of linear equations. 
Moreover, ( |36| ) could be expanded about points other than the origin, resulting 
in approximants accurate about these points. Our approach is to demand that 
(/-linear approximants shall reduce to linear approximants with polynomials of 
order nt , as q approaches unity, which corresponds to an approximation of the 
perimeter generating function by rational functions. Thus such approximants 
are only likely to be good if the perimeter generating function is dominated by 
a pole. We therefore expand (^) about t = and q = 1, taking into account 
only terms up to a fixed order N t (n t ) in t. In order to obtain the desired limit, 
N t (n t ) has to be chosen large enough. We order the resulting terms by increasing 
powers in 1 — q and, for each power, by increasing powers in t. For given N t {n t ), 
we compute g-linear approximants with polynomials of order n t in t and of order 
n q = 0, 1, . . . For given N t (nt), the highest obtainable order n q depends on the 
number of equations deriving from (^6|). To fix the multiplicative constant, 
we require ao(0,0) = 1 + ai(0,0). Information about the scaling function is 
obtained by applying the method of dominant balance to the approximants, as 
described above. 

In order to test the method of q- approximants, we will apply the method to a 
number of exactly solvable g-linear polygon models, which are mainly isotropic 
versions of the models defined above. These models have a rational perimeter 
generating function and can be shown to obey a g-linear equation of higher 



<k>(t,q)f(t,q) 



ai{t,q)f(qt,q) +b(t,q), 



(36) 
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than first order. Therefore, q-linear approximants (of first order) should give 
the correct differential equation for the scaling function. This can be checked 
by computing the leading coefficients in the asymptotic expansion of the model, 
using the anisotropic functional equation directly. If the functional equation of 
the isotropic model is known, the scaling function can alternatively be obtained 
by applying the method of dominant balance. 

We will illustrate the method of q-approximants by deriving the scaling 
behaviour of the q-linear models, defined above in the isotropic case, which do 
not obey a q-linear functional equation of first order. 

5.2 Stacks 

The generating function f(t,q) of stacks with equal horizontal and vertical 
perimeter activity x = y = t satisfies the g-linear functional equation of or- 
der 2 

ao(t, q)f(t, q) = ai(t, q)f(qt, q) + a 2 (t, q)f(q 2 t, q) + b(t, q), (37) 
where the polynomials ao(t,q), ai(t,q), a,2(t,q) and b(t,q) are given by 

a (t,q) = (l-q 2 tf(l-qt) 3 , (38) 

ai (t,q) = t 2 (l + q)(l-q 2 t) 2 , 

a 2 {t,q) = -q 3 t A , 

b(t,q) = -qt 2 {l-q 2 t){qH 3 - q 3 t 2 + q 2 t + qt-l). 

We found this relation by computing q-linear approximants of second order to 
stacks. It is possible to interpret this relation combinatorially Q]. The perimeter 
generating function f(t, 1) diverges at t c = (3 — \fh)j2 with a simple pole. The 
method of dominant balance can be applied to compute the scaling function 
about t c . It is of the same form as the scaling function for Ferrers diagrams, 
which corresponds to the observation made for the semi-continuous models ]l4| . 
We used stacks in order to test the method of g-linear approximants. Their 
rational perimeter generating function is obtained by approximants of cubic 
order (n t =3). A scaling analysis of cubic approximants at q = 1 and t c — 
(3— v5) /2 yields the correct type of differential equation for the scaling functions 
for each approximant, with generally incorrect coefficients. The accuracy of 
approximation increases with the degree n q in (1 — q) of the polynomials. For 
iV t (3) > 6 and n q > 1, the approximants yield the correct scaling equation. This 
result is robust against increasing the order n t . 

5.3 Decorated Ferrers diagrams 

We introduce Ferrers diagrams with a decorated top layer: Consider a decoration 
consisting of exactly k black squares. Approximants indicate that the scaling 
function is the same as for (pure) Ferrers diagrams. We checked this for k = 
0,1,. ..,5. 
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Instead of considering more complicated decorations of the top layer, we will 
now consider Ferrers diagrams with a decorated bottom layer. Let us approxi- 
mate the model with a decoration consisting of exactly two black squares which 
may be placed anywhere in the bottom layer. The model has the generating 
function 



G(x,y,q) = Y / ^- ^' ?™Z , (39) 



n(n- 1) y(qx) 



n=2 



2 {qy,q)n' 



where (t; q) n = Ofc=o (1 — 1 k t) denotes the g-product. The perimeter generating 
function of the model is 

We consider the isotropic case where x = y = t. The critical point is t c = 
1/2. Since the perimeter generating function is rational with numerator and 
denominator polynomials of order 4 and 3, we use q- linear approximants with 
m = 4 and Nt — 20. For n q > 2, the dominant terms of the approximants give 
the differential equation 

f + ±s) pfg) -(— + —) -P(s) - — = (41) 
16 J w \16 256/ ds y ' 128 v ' 

with exponents 

•-f, <t>=\- (42) 

This result is robust against varying the values of n q and nt in the approxima- 
tion. We checked this up to n q — 8 and also for increasing values of n q at nt = 5 
and n t = 6. The above equation leads to the scaling function 



''(*> = I / TT^yl dt ) 2e 85 (l + 16,- 2 ). (43) 

The asymptotic behaviour is given by 

2" 7 

P(s) ~ — (s^+oo) (44) 
P(a) - 2P(0)e 8s V (s-»-oo) 

no, - ^. 

Using the methods described at the end of the appendix, it can be shown that 
the isotropic model satisfies a g-linear equation of third order which can be used 
to test that the scaling function obtained by q- linear approximants is correct. 
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5.4 Ferrers diagrams with a 1-hole 

Ferrers diagrams with a 1-hole arc defined in the Appendix. We again consider 
the isotropic model where t = x = y. Since the perimeter generating function 
is rational with numerator and denominator polynomials of order 6 and 5, we 
use g-linear approximants with n t — 6 and N t — 20. For n q > 2, the dominant 
terms of the approximants give the differential equation 

( s 3 + —s) P(s) -(— + —) —PCS) - — = (45) 
V 16 J w \16 256 J ds K ' 128 v ' 

with exponents 

9 = 1, *=\ (46) 

This is the same equation as that for decorated Ferrers diagrams. This result 
is robust against varying the values of n q and rit in the approximation. We 
checked this up to n q — 8 and also by increasing values of n q at n t = 7 and 
nt = 8. 



6 Comparison with partial differential approxi- 
mants 

There is an existing approximating method designed to compute critical expo- 
nents and scaling functions about multi-singular points, known as the method 
of partial differential approximants (p.d.a.), due to Fisher and co-workers |7], [I?], 
|l5| , plj| |. We will briefly explain the method and compare it to our approach. 
The basic idea of p.d.a. derives from the observation that a scaling function 

,lt " = CT'(F?) + '" (47) 

obeys the first-order partial differential equation 

0f(t, q) + fo = 4>(t c ' t)^f(t, q) + (q c - q)^-f(t, q). (48) 

Therefore, p.d.a. of the form 

a(t, q)f(t, q) + b(t, q) = c(t, ?)^/(*. S) + ?)^-/(*. «), (49) 

at aq 

where a(i, g), b(t, q), c(t, q) and d(t, q) are polynomials in t and q, may serve to 
detect possible scaling behaviour about multi- critical points (t c , q c ) defined by 
the simultaneous vanishing of c and d, 

c(t c ,q c ) = d(t c ,q c )=0. (50) 

The critical exponents 9 and <f> can then be read off as lowest order coefficients in 
the expansion of the approximating polynomials about the multi-critical point. 
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Numerical methods can be used to determine subsequent terms in the Taylor- 
expansion of the scaling function. It has been shown Jl5[ ] that numerical inte- 
gration works if the crossover exponent cj> is restricted to 1/2 < cj) < 2. 

In contrast to this general setup, g-linear approximants are only suited to 
the detection of critical behaviour for models whose scaling function obeys a 
linear differential or difference equation of first order. Equivalently, g-linear 
approximants can be used to test whether a scaling function obeys an equation 
of the above type: If it does not, the approximants are likely to fail to converge. 
If it does, the approximants will converge and give the underlying differential or 
difference equation. This is then more specific information than can be gained 
from the p.d.a. approach, though the p.d.a. approach currently approximates a 
broader range of scaling behaviour. In subsequent work we will extend to non- 
linear approximants, which should combine the generality of the p.d.a. approach 
with the specificity of information obtainable by the g-approximants. 

We conclude with remarks about Ferrers diagrams and rectangles. Since 
Ferrers diagrams have a crossover exponent = 1/2, numerical integration 
using p.d.a. jL5| to obtain the coefficients of the scaling function is unlikely 
to converge. The model of rectangles (Q) does not obey a scaling law of type 
(p7|). It can be shown that the scaling function gives quadratic prefactors for the 
derivatives in (f48|). Therefore, the p.d.a. method cannot be successfully applied 
here in its standard form. For rectangles with the random chain considered 
above, p.d.a. may provide good estimates. 

In summary, we would expect that the method reported here is likely to 
be better than p.d.a. for those systems whose scaling function is described 
by, or well approximated by, a linear differential or difference equation of first 
order, while the p.d.a. method might be expected to better approximate those 
systems that do not. As we extend the method reported here to higher order 
g-functional equations, and possibly other types of functional equation as well, 
we would expect these new g-approximants to more appropriately represent a 
correspondingly larger class of systems. We emphasise that these are remarks 
of a general nature, and not the result of rigorous numerical comparisons, which 
we have not carried out. 



7 Conclusion 

We have developed techniques to obtain scaling functions for g-linear polygon 
models about the point where the perimeter generating function diverges, using 
the method of dominant balance. This led to scaling functions for a number of 
^-linear polygons models generalizing rectangles, Ferrers diagrams and stacks. 
The question arises as to what extent can the scaling behaviour for these models 
be obtained by direct methods? The most direct approach is to approximate 
the perimeter and area generating function by their Euler-Maclaurin sum and 
to estimate the resulting integral by methods of uniform asymptotic expansions 
[ fl9| , in the spirit of ]l^ . The authors are, however, not aware of standard 
methods to do this, apart from Bleistein's method which can be used to 
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analyze Ferrers diagrams. 

We introduced the method of g-lincar approximants to obtain scaling func- 
tions of models where an exact g-linear functional equation of first order does 
not exist. The method yields the correct scaling functions for exactly solvable 
models which can be described by g-linear functional equations of higher than 
first order, such as for decorated g-linear models. We claim that g-approximants 
will be appropriate for the analysis of the scaling behaviour of statistical mod- 
els whose scaling function may be well described by a difference or differential 
equation of first order. 

The method of dominant balance can also be applied to obtain differen- 
tial equations for scaling functions from g-functional equations different from 
g-linear. For example, it is possible to derive the differential equation for the 
scaling function of the model of staircase polygons which belongs to the 
g-quadratic class. This indicates that the idea of g-linear approximants can be 
generalized to more complex classes of polygon models satisfying g-algebraic 
functional equations. This leads to q-algebraic approximants. Even the q- 
quadratic class is interesting to analyze, since there are solvable models where 
the scaling behaviour is not known (such as staircase polygons with a hole; see 
also H for an example), and more interestingly it can be used to approximate 
the generating function for self-avoiding polygons, which displays a square-root 
divergence in its perimeter generating function. We are developing our method 
in that direction. 
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Appendix: Ferrers diagrams with a 1-hole 

We define the polygon model of Ferrers diagrams with a 1-hole and give a 
closed form for the perimeter generating function and for the perimeter and 
area generating function. 

The model of Ferrers diagrams with a 1-hole consists of all Ferrers diagrams 
where a unit square is removed from the interior. For a typical Ferrers diagram, 
this is graphically depicted in Fig. ||. If we denote the last column height by 
m, the last-but-one column height by mi, the bottom row length by n and the 
second row length by m, the number of possible sites for a hole h(n, m, m, mi) 
is given by 

h(n, m, ni, mi) = area — n — m — nx — mi + 4. (51) 
This translates into a formula for the generating function G(x,y,q) of Ferrers 
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Figure 3: A typical Ferrers diagram. A hole may occur anywhere in the shaded region. 



diagrams with a 1-holc 



G ( x > y> i) = - — — — ( F - < d - F ) - y( d v F ) + ^ F )) fa, qv, q), 



xy 

(1 - qx)(l - qy) 

(52) 

where F(x,y,q) is the generating function of Ferrers diagrams (without holes) 
of height and width greater than one. In particular, we have 

G(x,y,l)=( XV ) r. (53) 

\l-x-yj (l-x)(l-y) 

The perimeter generating function is rational with the same critical point as the 
Ferrers diagrams (without holes). The perimeter and area generating function 
is given by 

qx 2 y " tfyr ^ (k-2)q k x 

U [X ,y,q)- {1 _ qx)(1 _ qy) ^ iq2x . q)n 2^ x _ qhx , 154J 

where (t;q) n — YlkZoO- — q k t) denotes the g-product. G(x,y,q) satisfies a 
g-linear equation of third order. This equation can be derived by expressing 
G(x,y,q), G(qx,qy,q) and G(q 2 x,q 2 y,q) in terms of F, (d x F), (d y F), (d q F) at 

„3„. „3„ 



argument (q x, q y, q). This is done by using (|>2|) and the symmetrized version 
of the functional equation for F(x,y,q). The resulting system of linear equations 
can be solved for F in terms of G. Insertion of the result into the functional 
equation for F gives the functional equation for G. 
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